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ABSTRACT 

With several detections, the technique of gravitational microlensing has proven useful for 
studying planets that orbit stars at Galactic distances, and it can even be applied to detect 
planets in neighbouring galaxies. So far, planet detections by microlensing have been consid- 
ered to result from a change in the bending of light and the resulting magnification caused by 
a planet around the foreground lens star However, in complete analogy to the annual parallax 
effect caused by the revolution of the Earth around the Sun, the motion of the source star 
around the common barycentre with an orbiting planet can also lead to observable deviations 
in microlensing light curves that can provide evidence for the unseen companion. We discuss 
this effect in some detail and study the prospects of microlensing observations for revealing 
planets through this alternative detection channel. Given that small distances between lens 
and source star are favoured, and that the effect becomes nearly independent of the source 
distance, planets would remain detectable even if their host star is located outside the Milky 
Way with a sufficiently good photometry (exceeding present-day technology) being possible. 
From synthetic light curves arising from a Monte-Carlo simulation, we find that the chances 
for such detections are not overwhelming and appear practically limited to the most massive 
planets (at least with current observational set-ups), but they are large enough for leaving the 
possibility that one or the other signal has already been observed. However, it may remain 
undetermined whether the planet actually orbits the source star or rather the lens star, which 
leaves us with an ambiguity not only with respect to its location, but also to its properties. 
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1 INTRODUCTION 

Gravitational microlensing, i.e. the transient brightening of an ob- 
served star due to the bending of light caused by the gravitational 
field of an intervening foreground 'lens' star, was considered b y 
Einstein as early as 1912, as pointed out by iRenn et alj jl979l) . 
but he conclud ed that "there is no great chance of observing this 
phenomenon" ( lEinsteinlll936h . Only several decades of advance in 
technolo gy enabled the fi rst reported discovery of a mi crolensing 
event tAlcock et alll 19931) . following the suggestion by IPaczynskil 
( Il986h to use the technique as a tool for detecting compact matter 
in the Galactic halo. 

However, microlensing provides a valuable tool for a vari- 
ety of other astrophysical applications, and the most spectacular 
one nowadays is th e detection of extra-solar planets. It was al- 
ready pointed out bv lLiebe^ jl964) that the primary effect of plan- 
ets as gravitational lenses would be to produce a slight perturba- 
tion of the lens action of their respective host star. The distor- 
tion of the magnification pattern of the foreground lens star by 
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an orbiting planet, and the additional short blip or dip to the oth- 
er wise symmetric mi c rolens ing light curve was then first studied 
bv lMao &Paczvnskil jl99lh . A super-Jupiter was the first planet 
detected by this technique jBond et al] |2004|) . but its sensitivity 
even reaches below the mass of Earth, even for ground-based ob- 
servations (Bennett & Rhie 1996; Dominiketal. 2007). In fact, 
the detectability of planets below 10 has been impressively 
demonstrated with th e first discovery of a cool rocky/icy exoplanet 
teeaulieu et al.l2006l) . Microlensing is already singled out amongst 
all ground-based current techniques aiming at the detection of 
extra-solar planets by the respective host stars being at Galactic 
distances, rather than in the solar neighbourhood, and belonging 
to either of two stellar populations. Significantly beyond this, even 
planets orbiti ng stars in neighbouring galaxies, such a s M31, could 
be detected jC ovone et aljl2000l : Ichung et al.l l2006^■ whereas no 
other technique so far has been suggested that could achieve such a 
goal within foreseeable time. 

By creating a link between received light, the gravitational 
field of intervening objects, and relative transverse motions be- 
tween source, lens, and observer, the effect of gravitational mi- 
crolensing shows a substantial versatility. It is therefore not that 
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surprising that it provides other channels for detec ting planets or- 
biting stars other than the Sun. As suggested by I Graff & Gaudil 
ilOOd ). the light of close-in gas-giant planets would be detectable 
with large telescopes if such observations are scheduled while the 
planet follows its host star in exiting a caustic produced by a binary 
lens system, given that the light received from the planet would be 
far more strongly magnified than that received from its host star. 

Here, we discuss a further channel for revealing the existence 
of extra-solar planets from the study of microlensing light curves. 
Rather than considering a planet around the lens star, we study the 
effects of the orbital motion of the source star and a planetary com- 
panion around the common barycentre. While the planet is not seen 
itself, the small motion of its observed host star periodically al- 
ters the line-of-sight and thereby the relative lens-source position, 
which results in a change of the observed source magnification. In 
fact, an analogous effect is caused by the revolution of the Earth 
around the Sun, where the line-of-sight is altered due to the motion 
of the observer rather than the observed object, so that it constitutes 
a form of parallax effect. While t he annual para l lax in a microlens- 
ing event was first observed by lAlcock et alj ( 1 19951) . the orbital 
motion of stella r source b inaries has been studied extensively by 
vario us authors dOriest & Hu 1992; Han & Gould 1997; Paczyhski" 
Il997h . Some authors refer to the latter as 'xallarap' effect, but in 
fact, this nomenclature involves a double inversion, because par- 
allax is known as the apparent change in position of an observed 
object that is just the reflection of a change in position by the ob- 
server. Not surprisingly, both effects are not easy to distinguish, and 
a systematic analysis o f 22 microlensing parallax candidate events 
jPoindexter et alj|2005h explicitly finds 23 per cent of them being 
strongly affected by 'xallarap'. Future space missions capable of 
accurate astrometry such as SIM and GAIA will however allow 
measure the parallax along with the relative proper motion between 
lens and source star, allowing to obtain more accurate informatio n 
about the Galaxy (e.g. |Paczviiskill993 ; lRahvar & Ghassemil2005l) . 

We start our discussion in Sect. |2] with an introduction to 
the basics of gravitational microlensing and Keplerian motion, fol- 
lowed by a review of the annual Earth-Sun parallax based on the 
common form alism for orbital motion in observer, lens, and source 
developed by ( lDominil3ll998h . The discussion of annual parallax 
is subsequently transferred to the case of stellar reflex motion due 
to an orbiting planet, where the relevant parameters are identified. 
In Sect. [3] we turn our attention to the prospects for planet de- 
tection through our proposed channel, where we first discuss the 
strength of the observable effect and identify the favourable sce- 
narios. Subsequently, we discuss the input parameters and the re- 
sults of a Monte-Carlo simulation of synthetic light curves, thereby 
quantifying the planet detectability as a function of the planet mass 
and orbital parameters. Moreover, we study how to extract proper- 
ties of the planet and its orbit from the observed light curve. A short 
summary and final conclusions are presented in Sect.|4] 



2 PLANETARY ORBITS AND MICROLENSING 
2.1 Microlensing events 

According to the theory of General Relativity, a light ray passing 
beside a massive object is bent, the arrival time of its photons is 
delayed, and a ray bundle is distorted. Allowing for several possi- 
ble light rays from an observed source object to the observer, this 
phenomenon is commonly known as 'gravitational lensing'. If a 
foreground star with mass M at distance Dl happens to be suffi- 
ciently aligned with an observed background star at distance Dg, 



the time delay is negligible, while the angular separation between 
the images is of the order of micro-arcseconds, and therefore unde- 
tectable with current telescopes. However, contrary to extragalactic 
scenarios of gravitational lensing, stars in our own or neighbour- 
ing galaxies show a substantial proper motion, so that the image 
distortion results in an observable transient brightening as the fore- 
ground 'lens' star passes along the line-of-sight to the observed 
background source star, which is commonly referred to as a 'mi- 
crolensing event'. 

Explicitly, for lens and source being separated on the sky by a 
position angle u 9e, where the angular Einstein radius 



(1) 



provides the unique characteristic scale of gravitational microlens- 
ing, the combined magnification of the images reads jEinsteinl 
119M) 



A(u) = 



+2 



(2) 



uVu^ + 4 
where u — \u\. 

Let us consider source and lens star being in uniform proper 
motion with /Xg and /i^, respectively, so that the relative proper 
motion reads ^ — — /Xl ~ /i (cos tp, sin i/;). With uo marking 
the closest angular approach between lens and source star, realized 
at epoch to, the relative source-lens trajectory takes the form 



— sin ip 
cosi/j 



+ tJE — to) 



cos l/i 
sini/; 



where uje = ij./9e, so that 

ULs{t) = ^yu'^+UJ^{t~ toY . 



(3) 



(4) 



In the literature, rather than uj-e., a time-scale t^ = 

more frequently. Given that A{u) is strictly decreasing with u, 

and u{t) assumes a minimum of no at epoch to, the magnification 

AJWt)] leads to a symmetric light curve, peaking at to tPaczvriskil 

IT986h . 

The symmetry of the light curve is retained if the finite angular 
radius 6'* of the source star is taken into account. With a source size 
parameter = 6i,/6E and the abbreviations 
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I Witt & Mad d 19941) found a uniformly bright star being magnified 
by 
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for u 7^ p*, where K{k), E{k), and n(n; k) denote the com- 
plete elliptic integrals of first, second, and third kind, respectively, 
whereas 



A(m;p* 
for u = 



2 , 1 + 

P* pi 



arcsm ■ 



(V) 



p; + i 

p*. The size of the source star significantly affects the 
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observed light curve for u < 2 p*, while A{u; pi,) reveals Eq. ([2J 
as ^ 0. 

For Fa, denoting the intrinsic flux of the magnified source star 
and Fb that of the background, unresolved from the observed tar- 
get, the microlensing light curve is given by the received flux 



F{t)^FsA[u{t)-p.,I] + FB. 



(8) 



where the finite-source magnification A[u{t),pi,] is completely 
characterized by (uq, fo, cje, p*, I), where I stands for the bright- 
ness profile function of the source star, for which we assume 1=1 
(i.e. a uniformly bright source) throughout. 



2.2 Keplerian motion 

If one neglects relativistic effects, a planet and its host star of 
masses rrip and m*, respectively, both are in elliptic orbits around 
their common centre of mass. In fact, the motion of their separation 
vector r{t) = r^(t) — ri,{t) can be understood as a fixed virtual 
body of total mass 

m = TTip + m* (9) 

being orbited by another virtual body of the reduced mass 

m-r, rrii, 

M = — , (10) 

m 

and with rem it) denoting the motion of the centre of mass, 
rp(t) = remit) + {m*/m)r{t) 

r^t) = rcm(t) - (mp/m)r(t) (11) 

The orbit r{t) = [x{t),y{t),0] is characterized by its major 
semi-axis a, its orbital period P, its eccentricity e, the orbital plane, 
spanned by x (along major axis) and y (along minor axis), as well 
as the orbital phase at a reference epoch. Kepler's third law provides 
a relation between the orbital period P and the major semi-axis a 
by means of the total mass m, where 



P = 2n 



Gm 



(12) 



With positive x in the direction of periastron and positive y from 
periastron towards the motion of the planet, one finds 

x{t) = a[cos^(t)—e] , 

y{t) = ay/l-e^ sm^{t) , (13) 

where the eccentric anomaly ^ G [0, 2 7t) is given by 

^{t) - e sinC(t) = n{t - tp) - in{t - tp)J , (14) 

with n = (2 7t)/P, and tp being an epoch of periastron. 

For small eccentricities, an expansion in the lowest order of 
the eccentricity e is a fair approximation, which reads 



x{t) = a p{t) cos S,{t) , 
y{t) = ap{t)sm^{t), 

with 

p{t) = 1 - £ cos [n{t - tp)] 

and 

= n{t - tp) + 2e sin [n{t - tp)] 



(15) 



(16) 



(17) 



2.3 Earth-Sun parallax 

For the Earth, nip ^ m^,, and with the Sun at rest, Vp ~ i'{t). 
Moreover, the orbital major semi-axis is given by a = 1 AU, the 
orbital frequency is given by = (27r)/(l yr), and the eccentric- 
ity e = 0.017 is small. The effect of the annual parallax, i.e. the 
revolutio n of the Earth ar ound the Sun has been discussed in some 
detail bv lDominikl ( Il998h . and our subsequent discussion is based 
on the results derived in that context. 

The (parallactic) shift of the position of the observer {5r)o 
perpendicular to the line-of-sight is equivalent to a virtual displace- 
ment of the observed source star {5r)s- For the corresponding shift 
in the dimensionless angular coordinate {5u)s = {Sr)s/{Ds 0e), 
one finds 



(Su) 



Ds - Dl (<5r) 



With the relative lens-source parallax 



1 



(St) 



1 



(18) 



(19) 



it is customary to define a microlensing parallax tte = t^ls/Oe ■ 

A standard coordinate system based on the Earth's orbital 
plane is given by the ecliptical coordinates (/3, A), where positive 
f3 point towards ecliptical North, while increasing A follow the ap- 
parent motion of the Sun, which is in the same sense as the actual 
motion of the Earth. Moreover, A = is defined as the position of 
the Sun at vernal equinox. In fact, the ecliptical latitude /? equals 
the inclination of the Earth's orbit with respect to the line-of-sight 
to a source star at ecliptical coordinates A). 

In order to determine the effect of the orbital motion of the 
Earth, one needs to find its components perpendicular to the line- 
of-sight. While (x, y) span the orbital plane, where the positive 
X-axis points towards perihelion, the complementary z-coordinate 
points towards ecliptical North, thereby forming a right-handed 
three-dimensional system. We can also define a coordinate system 
with (ui, W2) in the plane perpendicular to the line-of-sight and U3 
pointing towards the observer, so that again a right-handed system 
is formed. Rather than using A as the longitude, it is more straight- 
forward to choose an angle if, where ip = corresponds to the 
Earth being at the perihelion. In fact, <^ = A + 7t + 1^3^, where if-, 
is the longitude of the vernal equinox as measured from the perihe- 
lion. While we can define ui (x z,U2 (x y, and 113 oc —x for /3 = 
and ((9 = 0, the coordinates for arbitrary {(3, ip) arise from a rota- 
tion around the fii-axis by the angle ip, and a subsequent rotation 
around the M2-axis by the angle — /3, so that 



fii(t) 

U2{t) 



1 AU 



7^(/3,<^) 



x(t) 

y{t) 



where 



n{f3,ip) = 



' sin f3 cos p — sin P sin ip 
— sin p cos p 



(20) 



(21) 



Since, in the small-eccentricity limit, Eq. l |15t , a rotation by p is 
equivalent to a shift in ^(t), one finds in general 



ui{t) = — TTE p(t) sin/3 cos ^(t) , 
U2{t) = TTE p{t) sin ({t) , 



(22) 



where ("(t) = ^{t) — p, and ^(t) being defined by Eq. (TT}. 

If we adopt the orientation angle of the source trajectory, 
given by iiLs(i). Eq. as referring to the (mi,M2) coordinate 
axes, we find for the total motion u{t) = WLs(i) + u{t), so that 
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its absolute square reads 

[u{t)f = + iJ^ {t - to)'^ + 

+ 2TT-E p{t) {sm({t) [uo cosi/; + uje (t — to) sini/;] + 
+ sin (3 cos ((t) [uo sin il^ — lje {t — to) cos V']} + 
+ ttI [p(t)f (sin' C(i) + sin' /3 cos' C{t)) ■ (23) 

With tiie Eartli's orbit defined with respect to the source at ecliptical 
coordinates (/3, A), only the microlensing parallax tte > 0, deter- 
mining the strength of the parallax effect, and the angle ^p, defining 
the orientation of the source trajectory, are free parameters along 
with {uo,to,u;E, p*) that define the magnification corresponding 
to the ordinary light curve including finite-source effects. Since, 
with this parametrization, uo refers to the minimal impact of the 
heliocentric trajectory, to no longer coincides with the maximum 
magnification. 

2.4 Stellar reflex motion due to orbiting planet 

Given that a planet and its host star orbit their common barycen- 
tre as described by Eq. Jilt , and that the motion of the observer 
and the source star are equivalent, it is obvious that the periodic 
displacement of the observed source star due to the orbiting planet 
and the annual parallax due to the revolution of the Earth around 
the Sun take exactly the same form. 

In analogy to the discussion of the previous subsection, let us 
define a parameter 



m Ds Oe 



>0, 



(24) 



which takes over the role held by tte in measuring the strength of 
the effect on the microlensing light curve. Let i denote the inclina- 
tion of the orbit with respect to a plane perpendicular to the line- 
of-sight, and let ip denote a longitude in the orbital plane that de- 
creases with the motion, where = at periastron. The position of 
the source star is then given by the orbital-plane coordinates {i, if), 
and similar to before, one finds with {5u)s ~ {Sr)s/{Ds Oe) 



ui{t) 

U2{t) 



XE T^/. X / x{t) 

— I yit) 



(25) 



with TZ given by Eq. | |2U . where the additional sign results from 
considering the motion of the source star rather than its planet. 

More explicitly, for small eccentricity e, one obtains with 
Eq. idD 



= XE p(i)sini cos^(i) , 
U2{t) = -XEp{t)smC{t), 



(26) 



with ^(t) = ^(t) — tf and ^(t) being defined by Eq. (TT}. Let us 
choose the coordinate axes of the motion of the barycentre WLs(i), 
as given by Eq. ([3}, as those of the (ui , M2 ) coordinate system. This 
yields a total motion u{t) = WLs(i) + whose absolute square 
is given by 

[u{t)f = ul+UJEit-to)^ + 

— 2 xe p{t) {sin ((t) [uo cos tp + LUE{t~ to) sin tp] + 

+ sin i cos (^(t) [uo sin ^ — oje — ^o) cos i/']} + 
+ xl[p{t)f (sin' C it) + sin' i cos' C (t) ) . (27) 

One indeed realizes that Eqs. i23\ and J27b have an identical form, 
where only tte <-> — Xe and (5 i. Correspondingly, uo refers 
to the minimal impact of the barycentre, so that again the epoch to 
does in general not mark a magnification maximum. 



However, in contrast to the effect of annual parallax, in ad- 
dition to the strength parameter xe and the direction of source 
(centre-of-mass) motion, characterized by the angle ip, there are 
a further 5 free parameters, namely the orbital eccentricity e, the 
orbital frequency fi, the orbital inclination i, the orbital longitude 
ip, and the time of periastron tp. In principle, the detection of reflex 
motion of the source star can be confused with the motion of the 
Earth around the Sun (the sign between tte and xe can simply be 
accounted for by + n), and we should take extreme care in 

avoiding to re-detect the habitable planet that we ourselves live on. 
It therefore needs to be ensured that at least one of the parameters 
determining the orbital motion is found to be incompatible with 
annual parallax. 

For circular orbits, i.e. e = 0, the time of periastron tp be- 
comes arbitrary, and we can choose to as reference epoch instead, 
so that along with pit) = 1, one can define C,{t) = ^1 (t — to) — f. 
In this case, besides (uq, to, 'i'E, p*), which describe the ordinary 
symmetric finite-source light curve in the absence of planets, the to- 
tal number of free parameters reduces to 5, namely (xe, 'ip,^l,i,(f). 

Another special case arises if the orbital plane is perpendic- 
ular to the line-of-sight, i.e. i = If this happens, modifying 
the angles iy3 or ?/; by the same amount has an identical effect, so 
that only one of them can be included in a set of independent free 
parameters, whereas the other is obsolete. 



3 PROSPECTS FOR PLANET DETECTION 

3.1 Strength of effect and favoured scenarios 

For a rough assessment of the amplitude of perturbations caused 
by the orbital source motion, and for an identification of the 
favoured scenarios, as well as on dependencies on system param- 
eters, let us first have a look at the strength parameter xe ~ 
irUp/M) [a/iDs Oe)], introduced in Eq. l l24t . At the source dis- 
tance Ds , the angular Einstein radius Oe corresponds to a physical 



Ds Oe 



= 2.8 AU 



'-^iDs-D.)^^ 



M 

TWr. 



1/2 



\Ds) 



(28) 



jQ / \ 1 kpc 

In accordance to what was found bv lHamadache et al] ( l2006h and 
ISumi et alj ( |2006|) . the strongest effects are therefore expected for 
small Ds — Dl, i.e. the lens stars being close to the source stars, 
in contrast to deviations by annual parallax, which are the strongest 
for lens stars close to the observer. This means that detections on 
stars in the Galactic bulge are dominated by lensing events caused 
by bulge stars rather than disk stars. A small Ds — Dl also implies 
Dl/Ds ~ 1, so that Ds Oe practically depends on the distance 
difference Ds — Dl only. Kepler's third law, Eq. l ll2t allows us to 
eliminate the major semi-axis a in favour of the orbital period P, 
so that 

-2/3 / \ -1/2 



XE 



6.4 X 10 



-4 "Ip 
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1 yr 



2/3 



Ds-Dr^-^/^ 



M 



0.3 Ma 



© 



1 kpc 



1/2 



(29) 



Please note that M is the mass of the lens star, whereas m = 
m-i, + nip is the total mass of the source system. While microlens- 
ing events prefer Af ~ 0.3 Mq, solar-mass source stars are much 
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brighter, so that those are the more reasonable target despite the fact 
that the mass ratio m^/m for a given planet mass is smaller than 
for low-mass stars. Moreover, solar-mass stars are far more likely 
to host gas-giant planets than low-mass stars. 
For the event time-scale ts, one finds 



tY. = 17d 



U60kms-iy 



M 



0.3 Mq 
1 



1 kpc 



1/2 



(30) 



while for Galactic bulge-bulge lensing, _Dl ^ ~ 160 kms~^ and 
Ds — -Dl ~ 1 kpc. Therefore, a characteristic value for the 
favourite scenario is given by ~ 17 d. 



3.2 Planetary signal 

Since displacements of the source star cause larger changes to its 
magnification the smaller the source-lens separation (and there- 
fore, the larger the magnification), the effects of the orbiting planet 
on the light curve increase with the source magnification, and the 
planet needs to be identified while the latter is substantial. There- 
fore, the characteristics of the planetary signal depend on whether 
the orbital period P is smaller or larger than the event time-scale Ie, 
as illustrated in Figs.[T]and|2] As by the definition of the parameters, 
uo does not refer to the closest approach between lens and source 
star, so that the peak magnification differs from Aq — u^^ — 125, 
where the discrepancy is larger for the long-period case with the 
larger xe. Given that we only altered the orbital period between 
the two cases shown, while leaving all other physical properties 
unchanged, the signal strength parameter xe increases with the or- 
bital period P, according to Eqs. ( 124b and l |29l ). 

For the case P <^ t^, the periodic pull of the planet on its 
host star leads to detectable ripples on the observed light curve, 
while the long orbital period P deprives us of such a characteristic 
signature for P ^ ts- Without a good indication of P, a distinc- 
tion with the effect of annual parallax becomes difficul t, and to the 
lowes t order, one only observes an acceleration effect iSmith et"ail 
l2003h . However, regardless of P/t-E, a best-fitting ordinary light 
curve leaves us with a mismatch near the tip of the light curve that 
is not overcome by adopting a different finite-source parameter p*, 
and allows a detection if the impact parameter uo is sufficiently 
small for such a signal to be prominent enough. 



3.3 Parameters of simulation 

After having identified the basic scenarios, we carried out Monte- 
Carlo simulations in order to study the detectability of planets or- 
biting Galactic bulge stars as a function of various parameters that 
describe the lens star, the source star, and its orbiting planet. For 
simplicity, we have assumed circular orbits (e = 0), so that the 
source magnification A[u{t); p*] is described by the 9 parameters 
{uo,to,LJE, p*, XE, i/', i, <p)- 

Without loss of generality, we set to ~ 0. Moreover, we con- 
sider planets of mass nip — 1 Afjup or nip = 10 Mjup orbiting 
a source star of m* = 1 Mq. With these choices, according to 
Eq. i29h the strength parameter xe then becomes a function of the 
lens and source distances Dl and Ds, for Ds — -Dl ^ Dl essen- 
tially of Ds — Dl, the orbital frequency Q — {2n)/P, where 
P denotes the orbital period, and the lens mass M. In relation 
to the sampling interval and the event time-scales, we generate a 
uniform distribution in lg[P/(l d)] ranging between 3 d and 1 yr. 



Source star microlensing wobble 

caused by orbiting planet 

mp = 10 Mjup, m^, = 1 Mq, P = 50 d 
tE = 80 d, MO = 0.008, p* = 0.0036 
Xe = 2.0 X IQ-^, i = 85°, = 280°, ip = 
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Figure 1. A model light curve for which the orbital period P is smaller 
than the event time-scale = ~ ^ . The adopted parameters are compat- 
ible with a planet of mass mp = 10 Mjup orbiting a star of solar mass 
and solar radius at distance Ds = 8.5 kpc, whereas the lens stai' of mass 
M = 0.4 Mq is located at Dl = 8.0 kpc from the observer The upper 
panel shows the resulting light curve (solid) as well as a best-fitting approx- 
imation (dashed) with an ordinary model {uq, to , u>e , p*) ■ The difference 
between these two curves is plotted in the lower panel. Dotted lines indicate 
deviations by 2 per cent or 0.3 per cent, respectively. 



While we adopt a 'natural' uniform distribution of impact param- 
eters Uo G [0, 1], not taking into account any selection bias by the 
experiment, the phase angle ip, orientation angle tp, and inclination 
angle i are all naturally uniformly distributed, where ip,ip & [0, 2tt) 
andi G [0,7r/2]. 

The event time-scale Ie follows from drawing a lens dis- 
tance Dl, source distance Ds, velocity v, and lens mass AI from 
the adopted distributions for the Galactic bulge described in Ap- 

'q. 



pendix |A] Moreover, with assuming a source radius i?* = 1 P, 
one obtains the source size paramete r p<, = \Rt/(v tE) ] (Dl/Ds). 

In contrast to the discussion bv lRahvar et al.l ( l2003i) of the ob- 
servability of parallax effects towards the Large Magellanic Cloud, 
which used parameters specific to the EROS (Experience de la 
Recherche d'Objets Sombres) campaign, we adopt the simple prag- 
matic approach of assuming a constant photometric uncertainty for 
equally-spaced observations over the course of the microlensing 
event without any loss due to bad weather. In particular, we choose 
(a) a sampling interval of 2 hours with 2 per cent accuracy, re- 
sembling current follow-up observations, or (b) a sampling interval 
of 15 minutes with 0.3 per cent, resembling upcoming campaigns. 
More precisely, we demand the fractional uncertainty of the mag- 
nification to match the quoted value, or equivalently, the flux af- 
ter subtraction of the background to be measured that precisely. In 
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Figure 2. A model light curve for which the orbital period P is larger than 
the event time-scale t^. The adopted parameters are compatible with the 
same scenario as chosen for Fig. [T] As before, the upper panel show the 
resulting light curve (solid) together with a best-fitting ordinaiy light curve 
(dashed), while the lower panel displays their difference. Dotted lines refer 
to deviations by 2 per cent or 0.3 per cent, respectively. 



fact, observing campaigns need to account for such a requirement, 
or time is being wasted on taking data on strongly blended targets 
without a chance to extract meaningful results. The simple choices 
allow us to focus on the primary dependencies without being bound 
to the variety of different existing or possible setups, which all show 
different characteristics with respect to the crowding of targets, the 
distribution of event impact parameters for main-sequence stars, 
and the actually achieved photometry, where all these effects are 
correlated with each other. 



3.4 Detectability in simulated events 

For each of the created synthetic light curves, we obtain best- 
fitting model parameters by means of -minimization, which cor- 
responds to a maximum-likelihood estimate for Gaussian error 
bars, for an ordinary light curve (uq, to, ^^e) and independently the 
8 parameters («o, io, f^E, Xe, V*! *i ¥') that include the descrip- 
tion of the motion of the source induced by a planet in a circular 
orbit. In both cases, we also determine effective best-fitting source 
and background fluxes _Fs and Fb, taking into account a potential 
difference between the best-fitting magnification and the true mag- 
nification arising from the simulation. We explicitly adopt a finite 
angular radius 8^ of the source star, but do not refit for the respec- 
tive parameter p* = 9-^/9^, in order to save computing time and to 



avoid parameter degeneracies if finite-source effects are not promi- 
nent. 

A quantitative measure for the detectability of a planetary sig- 
nal then results from a likelihood-ratio test, involving the respec- 
tive minima, namely (Xmin)^°' and (Xmin)'''''"°'- In fact, with 
AXp = (xfnin)'"' - (Xmin)'"""^' following a "distribution with 
5 degrees of freedom, one finds a probability P(Axp ^ 11.07) = 
0.05 for such a difference to arise. At this significance level, we 
therefore reject the hypothesis that an ordinary light curve explains 
the data whenever Axp ^ 11.07, and claim the detection of a sig- 
nal. 

A substantial fraction of such 'detections' however do not in- 
volve characteristic features, so that instead of the assumed effect of 
periodic source motion revealing the presence of an orbiting planet, 
these could be of different origin if found in an observed event. In 
order to account for false positives, we therefore adopt a second 
criterion. In analogy to the primary criterion, we determine the 5 
parameters (ito, to, we, ttb, '(/') that include the annual parallax due 
to the Earth's revolution, and carry out a similar likelihood-ratio test 



based on the criterion Axg 



(y2 . -.(0) _ . ) 



^ 5.99, 



where the adopted value results from P(Ax© ^ 5.99) = 0.05 for 
a x^-distribution with 2 degrees of freedom. We then consider the 
planetary signal to be characteristic if the detection of a signal is 
significant (Axp ^ 11.07), while a signature of annual parallax 
is not (Ax® < 5.99); and the signal to be featureless otherwise. 
With this procedure, we not only reject events with deviations that 
are likely to be caused by the motion of the Earth, but in an ele- 
gant way, we also get rid of all uncharacteristic deviations that are 
compatible, which could e.g. be due to differences in the optimal 
finite-source parameter or weak effects resulting from binarity of 
the lens or source object. At the end, we are left only with cases 
for which either the observed deviation is characteristic for the pe- 
riodic motion of the source star, or where an asymmetric feature is 
closely mimicked by effects other than the annual parallax. 

The relative abundances of planet detections amongst all sim- 
ulated events as a function of various parameters, namely the event 
impact parameter ito, the orbital period P = {2n)/Q., the event 
time-scale = lo^^, the lens-source distance Ds~Di^, the orbital 
inclination i, and the resulting strength parameter xe are shown in 
Figs. [3] and |4] corresponding to planets of mass nip = 10 A/jup 
or nip — 1 Mjup, respectively, orbiting solar-mass Galactic bulge 
stars, and the two adopted observational capabilities with regard to 
the sampling frequency and achievable photometric accuracy. 

For nip = 10 Mjup and the less favourable 2 per cent ac- 
curacy and 2 h sampling, about 1/4 of the significant deviations 
are characteristic for a detection of orbital stellar motion against 
false positives, providing a detection efficiency of ~ 2.5 per cent 
over the adopted sample. Small impact parameters strongly sup- 
port the detection, with efficiencies reaching ~ 24 per cent for 
characteristic signals arising amongst the smaller number of events 
with ito < 0.02 (Ao > 50). While nearly all detected signals are 
uncharacteristic for uo > 0.1, characteristic and uncharacteristic 
detections are of similar frequency for smaller ito. The detection 
efficiency is found to increase towards larger event time-scale 
and orbital period P, as well as to favour smaller Ds — I5l and 
face-on orbits i ~ 90° over edge-on orbits i ~ 0° . 

With a close-to-continuous sampling at 15 min intervals and 
0.3 per cent photometric accuracy, planets of 10 A/jup become hard 
to miss, with more than 85 per cent of them showing detectable sig- 
nals, and a small uo no longer being in strong favour of a detection. 
As compared to the less favourable observational capabilities, more 
than half of the detected deviations contain characteristic features. 
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Figure 3. Efficiency of detecting a planet of mass mp = 10 Afjup by means of the effect of the orbital motion of its observed host star in the Galactic bulge 
with m* = 1 Mq on the light curve arising from its light being bent by an intervening foreground Galactic bulge lens star as a function of the strength 
parameter XE, defined by Eq. (24), the event impact parameter G [0, 1], the event time-scale tE, the orbital period P, the lens-source distance Da, — Dj^, 
and the orbital inclination i. For the left panel, we assume observations with 2 per cent photometric accuracy regularly spaced at 2 h intervals, while the results 
shown in the right panel correspond to increased observational capabilities of 0.3 per cent photometric accuracy and a sampling interval of 15 min. The hatched 
area refers to cases where the signal is characteristic enough to be distinguished from false positives. 




Figure 4. Detection efficiency for a Jupiter-mass planet (mp = 1 Afjup) from the influence of source-star orbital motion on the microlensing light curve for 
two different observational capabilities, similar to Fig.|3] but considering only events with smaller impact parameters mq G [0, 0.1]. 
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This is a substantial improvement, pushing the detectabihty of plan- 
ets by a factor of ~ 20. The much larger gain on the detectability 
with the improved capabilities for shorter periods accounts for a 
decrease of the characteristic detections towards larger orbital peri- 
ods (whereas the fraction of uncharacteristic detections increases). 
In particular, for all the shorter orbital periods P < 10 d, the plan- 
etary signal is characteristic, and almost 60 per cent of the planets 
are found. One might be puzzled by the fact that for the largest en- 
countered strength parameters xe, the prospects for a characteristic 
detection decrease. While it seems a bit surprising, this behaviour 
results from xe being correlated with the orbital period P, and 
larger orbital periods mean that the light curves contain less char- 
acteristic features. 

For the prospects of detecting a Jupiter-mass planet, the results 
of our simulations reveal mostly the same trends with the param- 
eters as for 10 A/jup. While for such less massive planets, one is 
extremely unlikely to succeed unless uo < 0.03 with a strategy of 
2 h sampling at 2 per cent photometric accuracy, and a substantial 
rise of the detection efficiency towards smaller uo even for the more 
favourable observational setup, very small uq < 0.0003 are not op- 
timal either, given that finite-source effects wash out signals in that 
regime. On average, characteristic detections happen for ~ 0.9 per 
cent of all events with uo < 0.1, but for ~ 12 per cent of events 
with Uo ^ 0.005. With the adopted better observational capabili- 
ties, these values rise to ~ 30 per cent or ~ 70 per cent, respectively. 

Given that the same value of the strength parameter xe for a 
larger planet mass nip implies a smaller orbital period P, which 
eases the detection of a planetary signal, the detection efficiency 
for more massive planets at same xe should be larger. However, 
the smaller detection efficiency for some of the ranges shown in 
the figures is a result of averaging over events with uo G [0, 0.1] 
for rUp = 1 Mjup, whereas an average over uq G [0, 1] has been 



taken for nip 



10 Mjup. 



Rather than the planet detection efficiencies for certain pa- 
rameter ranges. Figs, |5] and |6] show the distribution of the detec- 
tions with these parameters, i.e. in contrast to before, the distribu- 
tion of the respective parameter ranges among the simulated events 
has been considered (which makes no difference for uniformly dis- 
tributed quantities, which have been omitted). In all four consid- 
ered cases (two planet masses and two sets of observational capa- 
bilities), the fact that microlensing events preferentially arise with 
time-scales 5 d ^ tE ^ 80 d leaves such small with the largest 
number of planet detections, despite the fact the the detection effi- 
ciency increases with t^- With 2 h sampling and 2 per cent accu- 
racy, there are different preferences for the impact parameter uo, 
depending on the mass of the planet. While for rUp — 10 A/jup, the 
bulk of detections is expected from 0.01 ^ uo ^ 0.3, the difficulty 
of revealing a signal for larger uo in the case of for rUp = 1 Mjup 
moves the preferred range to 0.0003 < uo ^ 0.03. For both planet 
masses, the better observational capabilities of 15 min sampling 
and 0.3 per cent photometric accuracy lead to the larger number of 
events with larger uo providing the larger number of detections. 



3.5 The challenge of detecting extragalactic planets 

No other technique but gravitational microlensing has so far 
been suggested to achieve the detection of planets orbiting stars 
in neighbouring gal a xies, such as M31, within foreseeable time 
jCovone et al] I2OO0I : IChung et all l2006h . Given that the signal 
strength xe related to the periodic reflex motion of the source star 
due to an orbiting planet lacks of practical dependence on the dis- 



tance Ds, one might hope that this planet detection channel could 
also be a viable alternative for extragalactic planets. 

As compared to the Galactic bulge, there are however some 
severe difficulties that make any such attempts a quite demanding 
challenge. Most importantly, one faces the problem that the target 
stars are not resolved, so that the effectively achieved photomet- 
ric accuracy is quite limited. Unless major technology leaps are 
made, one is therefore restricted to very strong signals, which are 
unlikely to occur. Catching main-sequence source stars rather than 
giants will only be possible during phases of extreme magnifica- 
tion, limiting the number of suitable events further. Beyond that, 
detections on M31 targets are disfavoured due to the larger typi- 
cal Ds — D\, ^ 10 kpc as compared to the Galactic bulge, where 
Ds - -Dl ~ 1 kpc. 

3.6 Does the planet orbit the source star? 

So far, we have assumed by construction that the periodic alteration 
of the relative position between lens and source star as seen from 
Earth is the result of a planet orbiting the observed source star, and 
only the revolution of the Earth itself has been considered as an 
alternative. However, as explicitly pointed out to us by V. Bozza 
(private communication), a similar effect might also arise from a 
planet orbiting the foreground lens star rather than the observed 
background source star. 

We are quite familiar with planets orbiting the foreground lens 
star revealing their presence due to them affecti ng the gravitational 
bendi ng of light received from the source star jMao & PaczYhskj 
Il99lh . Notably, the arising characteristic signal represents a snap- 
shot of the planet at its current angular separation from its host star, 
and does not depend on its orbital period. Moreover, the chances of 
detecting a planetary signal rely on a resonance of the angular Ein- 
stein radius 9e with the angular separation of the planet 6p — d Se, 
whereas the planet is likely to escape detection for either d ^ 1 or 
d> 1. 

It is obvious that similar to a source star of mass m, a lens star 
of mass Ad also exhibits a periodic shift due to a planet orbiting at a 
semi-major axis a, whose strength in full analogy to the discussion 
of Sect. |2.4| is given by 



Ae = -TT 



Inserting characteristic values yields 

M 



(31) 



1.4 X 10" 



P 
1 yr 



2/3 



1 kpc 



0.3 Mq 

1/2 



-7/6 



-1/2 



(32) 



and by comparing this expression with the corresponding Eq. J29b . 
one sees that the detection as compared to planets orbiting the 
source star is not only facilitated by the smaller typical stellar mass, 
but lens distances not only close to the source star, but also close to 
the observer provide favourable configurations. 

Assuming circular orbits, ajypical orbital radius is given by 
a ~ 1.36 d Di, Oe, ( lDominikl2006h . so that for d — 1, characteristic 
orbital periods are of the order of 



Po = 5.6 yr 



Dl-Ds 
1 kpc 



3/4 



\Ds 



3/4 



M 



0.3 Mq 



1/4 



(33) 



which usually exceed those for which we expect to detect signa- 
tures of the period reflex motion. Given that Po oc d'*^^, we find 
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Figure 5. Distribution witii tlie signal strengtli parameter XE. the event impact parameter uo G [0, 1], the event time-scale tE, or the lens-source distance 
Dg — Di^ of the detection of a planet of mass nip = 10 Mj^p orbiting a Galactic bulge star with m* = 1 Mq from the change in the microlensing light 
curve resulting from the periodic displacement of the observed star caused by the gravitational pull of the orbiting planet. As for Fig. [3] the two panels refer to 
different observational capabilities, namely 2 h sampling with 2 per cent accuracy or 15 min sampling with 0.3 per cent accuracy, and the hatched area marks 
the characteristic detections. The total area of the bins gives the average detection efficiency. For the uniformly-distributed quantities lg[P/ (1 yr)] and i, the 
distribution of the detections is proportional to the detection efficiencies plotted in Fig. [5] 
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Figure 6. Distribution of the planet detections with several parameters, similar to Fig.|3] but now for a planet mass of nip = 1 Afjup and the impact parameter 
in the range ug G [0, 0.1]. 



in particular that orbital periods Po ^ 50 days roughly correspond 
to separation parameters d < 0.08, so that substantial deviations 
due to gravitational bending of light are not likely to occur. In fact, 
such can always be checked for by means of constructing model 
light curves that take this effect into account. 



However, as a consequence, we are facing the situation that 
even in the case that one identifies a characteristic signal that cannot 
arise from the revolution of the Earth, it remains possible that one 
ends up with two competing interpretations putting the planet in 
orbit around the source star or the lens star, respectively. 



3.7 Properties of the planet and its orbit 

For the 'usual' channel of planet detection by microlensing, the 
properties of the planet and its orbit affect the light curve only by 
means of two dimensionless parameters, which can be taken as the 
planet-to-star mass ratio q and the separation parameter d, where 
d 6e is the instantaneous angular separation of the planet from its 
host star. With the technique being most sensitive to planets around 
d ~ 1, and planetary signals lasting between hours and days, these 
are substantially smaller than the orbital period P, and therefore 
the only information about the planetary orbit is provided by the 
separation parameter d, whereas neither the orbital eccentricity e, 
or the orbital inclination i can be inferred from the snapshot. 

The mass of the planet Mp = g Af is related to the mass 
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of the lens star M, which does not follow directly from the light 
curve, but in general needs to be inferred probabilistically from 
the event time-scale tE ~ Oe/i-i, with the angular Einstein radius, 
given by Eq. (TJ, being a function of the lens mass as well as of 
the lens and sou rce distances Dl and Ds ■ As discussed in detail by 
iDominikI j2006l) , this requires the adoption of a kinematic model 
of the Galaxy as well as of mass functions of the underlying stellar 
populations that make up the lens stars. Light curves that involve 
planets with small masses however are likely to be influenced by 
the finite angular size of the source star. With the possibility to 
determine 0* from stellar typing based on its magnitude and colour, 
and better with a spectrum, the extractable time-scale ti, — 
in which the source moves by its own angular radius with respect to 
the lens, allows to infer the proper motion jj,, and thereby with the 
event time-scale tE ~ ^E//iof the angular Einstein radius 9e ■ With 
a reliable estimate for the source distance Ds, the mass of the lens 
star M thereby only becomes a function of the lens distance Dl, 
reducing the uncertainty substantially. However, the mass of the 
lens star M, and thereby the mass of the planet Mp = qM only 
follows from the observed light curve, if moreover the microlens- 
ing parallax tte ~ ttls /^e can be determined. As for the lens mass 
A/, only a probability density can in general be obtained for the 
instantaneous physical projected separation f = dDi^ 6e, while a 
stochastic distribution for the orbital semi-major axis a further fol- 
lows with an orbital deprojection and assumption of a distribution 
of orbital eccentricities e. Similarly, probabilistic estimates for the 
orbital period P = 2 7t [a^/iGM)^/'^ can be derived. 

In contrast to the just 2 additional parameters (d, q) as com- 
pared to an ordinary microlensing light curves, our alternative 
channel of detecting planets from the orbital motion of observed 
microlensing source stars around the common barycentre involves 
7 parameters (xe, i/), il, i, e, tp). For studies of the planet pop- 
ulation, the values of the inclination z, the phase angle (p, and the 
time of periastron tp are of little use. Moreover, while for the annual 
parallax, the direction angle of the source trajectory with respect 
to the ecliptic plane, thereby being well-defined in space, carries 
useful information, the direction angle with respect to the orbital 
plane of the planet does not. In sharp contrast, a direct determina- 
tion of the orbital period P = (2 7t) /SI and the orbital eccentricity 
e are valuable. 

Since spectral typing not only provides an estimate of the an- 
gular source size S*, but also of the stellar mass m*, the orbital 
semi-major axis a = Gm[P'^ /{in^)]^^^ is determined alongwith 
the orbital period P, given that m « m*. The mass of the planet 
therefore follows as nip — xe [{Ds dE)/a] m, where the main 
indeterminacy results from the unknown Oe with in general only 
tE = ^e/m being extractable. However, with finite-source effects 
being observed, 6e is measured, and thereby the planet mass nip 
will result. Given the fact that we do not require any knowledge 
about the lens distance Dl, a measurement of the microlens paral- 
lax tte ~ ttls /Oe is not helpful in this case. Otherwise, as for the 
'standard' channel, one only finds a probability density for mp, but 
in contrast to that case, one obtains measurements for the orbital 
semi-major axis a, the orbital period P, and the orbital eccentricity 
e. 

While the parameters (xE,tp, i, f, £, tp) are extractable in 
principle, the power for determining the properties of the planet 
and its orbit are limited by severe degeneracies that occur in sev- 
eral (not unlikely) cases. For example, a proper measurement of 
the orbital period P is only possible for P < StE, whereas for 
wider orbits, the effect on the microlensing light curve is mainly 
described by the acceleration of the source trajectory in the vicinity 



of the ep och at which the or dinary, unperturbed light curve reaches 
its peak jSmith et al.ll2003l) . Moreover, for tE < P < 3tE, the 
orbital period P appears to be strongly degenerate with the orbital 
inclination i and cannot be properly disentangled. 

Finally, with a large variety of effects causing small devia- 
tions, one needs to take care to investigate all possible alternatives 
such as lens or source binarity, and in particular the annual Earth- 
Sun parallax. As already pointed out in Sect. 12.41 the latter causes 
an identical signature, with the only difference that the parameter 
space is restricted to the fixed values that define the Earth's orbit. 



4 SUMMARY AND FINAL CONCLUSIONS 

We find that an alternative channel for the detection of extra-solar 
planets by microlensing is provided by the orbital motion of the 
source star around the common barycentre with an unseen plane- 
tary companion, as opposed to the standard channel where a planet 
orbiting the lens stars alters the bending angle and thereby the ob- 
served magnification. We derived a formalism in exact analogy to 
the treatment of the annual parallax that results fr om the revolutio n 
of the Earth around the Sun, mainly following iDominikI ( Il998h . 
which however involves 7 additional free parameters as compared 
to an ordinary microlensing light curve, which reduce to 5 for cir- 
cular orbits, whereas there are just 2 for the annual parallax due to 
the known properties of the Earth's orbit. 

Constituting an exchange of the roles of source star and ob- 
server as compared to parallax effects, which are most prominent 
for lens stars close to the observer, the strongest signals arise for 
lens stars close to the source star. In this limit, the signal strength 
practically depends on the difference Ds — Dl of the source and 
lens distances only, rather than on the individual values, so that 
planetary signals on observed Galactic bulge stars will show up 
predominantly in events that involve lens stars in the Galactic bulge 
rather than the Galactic disk. Moreover, the signal strength persists 
for source stars in neighbouring galaxies, such as M3 1 , but achiev- 
ing a sufficient photometric accuracy (on an unresolved) target for 
being able to claim a detection is extremely challenging. 

Other than for the standard microlensing channel, which is 
strongly biased in favour of K- and M-dwarfs due to their larger 
abundance, more massive stars are the most prominent targets due 
to their greater brightness, which results in favourable prospects for 
studying gas-giant planets, which are known to be extremely rare 
around low-mass stars. 

Unless the orbital period can be identified from the observed 
light curve, there are a variety of alternative explanations for the 
nature of the event that are compatible with the acquired data. 
Rather featureless deviations might also result from finite-source 
effects, the revolution of the Earth around the Sun, as well as lens 
or source binarity. Therefore, a proper characterization can only 
be expected if the orbital period P does not substantially exceed 
the event time-scale tE, which works against the fact that signals 
increase with P. Apart from this, planets in closer orbits, in partic- 
ular with P < 50 d, could be as well orbiting the lens star rather 
than the source star, given that their effect on the bending of light 
would be expected to be negligible. 

A rough estimate of the signal strength showed that the ob- 
servability of signals with current experimental setups is practically 
limited to massive gas giants, and a Monte-Carlo simulation of a 
survey with 2 h sampling and 2 per cent photometric uncertainties 
revealed that, for the Galactic bulge, the detection of Jupiter-mass 
planets will be dominated by events with small impact parameters 
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Mo ^ 0.03, whereas there is a substantial chance to detect plan- 
ets of mass TTip ~ 10 A/jup for the larger sample of events with 
fto ^ 0.1. Microlensing observations with improved capabilities, 
namely 15 min sampling with 0.3 per cent accuracy, would sub- 
stantially increase the prospects, where the detection efficiency for 
Jupiter-mass planets in events with uq < 0.1 being boosted from 
0.9 per cent to 30 per cent, whereas the prospect for detecting plan- 
ets that are 10 times more massive amongst all events with uo < 1 
is pushed from 2.5 per cent to 50 per cent. With the much higher de- 
tection efficiencies for moderate impact parameters, upper decades 
of Uo (which carry more events) would provide a larger fraction of 
the detections than lower decades. 

If one is able to properly measure the orbital period P from 
the microlensing light curve, which requires P < t^, the orbital 
semi-major axis a can also be determined provided that typing of 
the source star yields its mass m*. Along with the orbital eccentric- 
ity e, substantially more information can be extracted as compared 
to the standard channel that lacks of vital constraints on the plane- 
tary orbit. The planet mass nip however still depends on the angular 
Einstein radius 6e, which is not known in general, leaving the need 
to adopt a kinematic model for the lens and source populations, as 
well a s a mass function for the lens stars, as discussed by Dominik 
( l2006h . This can only be overcome by additional measurements of 
either 6e or the lens-source proper motion /i, e.g. from the assess- 
ment of finite-source effects involving the time-scale = 6**//^, 
where the angular size of the source star follows from typing 
alongside the mass m*. 

Current microlensing observing campaigns are not expected to 
provide many such detections, but with a detection probability of a 
few per cent for ni p = 10 Mjup, and an estimated abundance also 
of a few per cent l lUdrv & Santos! 12007^1, a corresponding signal 
may already be present among the several thousand events that have 
been monitored so far (with 700-1000 new ones currently being 
discovered every year). However, finding it may require a careful 
extensive data-mining effort, given that deviations of this kind are 
easily missed or misidentified. 
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APPENDIX A: MODELS OF THE GALACTIC BULGE 
POPULATION 

Given that the number of source stars in the observed field is pro- 
portional to ps{Ds) Dg, while the event rate for a given source star 
is proportional to pl(£'l) Dl Oe, where ps and pL denote the vol- 
ume mass densities of the source and lens stars, respectively, a joint 
probability density for the source and lens distances is proportional 
to the differential total event rate 

j2t 



d^r 



oc Pl(-Dl) ps(£>s) X 



xoi^i^L (i--p^) e(Ds-OL)e(DL) (Al) 

with Q{x) denoting the step function. This however does not ac- 
count for the fact that more distant source stars appear fainter on 
average and that their light is more likely to be affected by extinc- 
tion. Integration of Eq. JAIb over the source distance Ds leads to 
the probability density of the lens distance Dl > being propor- 
tional to 

dr 



Ds 



dDs. 



(A2) 



oc pLiDL)VD^ J ps(Ds)Dsa/1 

Alternatively, the joint probability density as given by Eq. iAll can 
be transformed to refer to the source-lens distance Dls = Ds — Di^ 
rather than the lens distance Dl, yielding 

,n Pl(I>s - I5ls)ps(I5s) X 



xDI 



Dls ( 1 - ^ ) ©(Dls) e{Ds - Dls) . (A3) 



Therefore, one finds the respective probability density of Dls as 

oo 

dr 



dDLs 



/Dls J Pl(Ds -Dls)ps(Ds) X 



1 ^LS , n 

1- — dDs. 



(A4) 
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Figure Al. Joint probability density V-^\A^V j (ADi, dDs)] (Dl,-Ds), 
indicated by equidistant contours between zero and the maximum, of the 
lens distance Dl and the source distance Ds for self-Iensing withing the 
Galactic bulge in the direction of i = 0° and h= —3°. 
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Figure A2. Probability density of the lens-source distance 
(dr/dDLs) (i'Ls) for bulge-bulge lensing, where the source 
stars have been located towards I = 0° and 6 = — 3° and the bulge mass 
density is given by Eq. (A5). 



Amongst the various laws suggested for the mass density of 
the Galactic bulge, we adopt a triaxial barred shape of the form 

A^bulgo 



Pbulgc 



where 



6.57TTSxSySz 



exp I 



V2) 



+ (-r 



(A5) 



(A6) 



which is a favourite for matching COBE data i lDweketal.lll995t) 
with Sx = 1.49 kpc, Sy — 0.58 kpc and = 0.40 kpc, and 
Mbuigo = 1.7 X 10^° Mq. While 2 measures the distance towards 
Galactic North, the axes referring to the x and y coordinates are 
tilted by an angle 9 ~ 20° with respect to the direction towards 
the Galactic centre and the direction of the local circular motion, 
respectively. 

We neglect the variation with the sky coordinates of the source 
star, and adopt a typical average Galactic longitude / = 0° and 
latitude b — —3°, so that 



X = {Dl cos b — Rcc) cos 8 
y = — (Dl cosfe — -Rgc) sinS 
2 = Dl sin b , 



(A7) 



where Rgc ~ 7.62 kpc denotes the distance to the Galactic centre. 

Given that the bulge self-lensing is dominant, pbuige describes 
both the lens and the source stars. The joint probability density 
of the lens distance Dl and source distance Ds, according to 
Eq. iAH . as well as the probability density of the lens-source dis- 
tance _Dls = -Ds — Dl, according to Eq. l lA4t , respectively, that 
result from this choice, are displayed in Figs. lAll and lA2l 

For the velocities of lens and source stars, we adopt a 
two-dimensional Maxwell-Boltzmann distribution in the absolute 
transversal velocity v, where the probability to find v in the inter- 
val [v, V + dv] is given by 



>l>„(i;) dv = — — ■ exp 
2 



V 



dv . 



(A8) 



While we adopted cr = 100 km s , we neglected its effective vari- 
ation with Dl and Ds , as well as the motion of the Sun, which are 
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Figure A3. The adopted normalized Galactic bulge mass function 
'^Is(m/Mq)[^s{^/^q)]' which corresponds to the probability density 
in lg(M/M0) for the mass of a given star. 



fair approximations for Dl/Ds ~ 1, marking the region dominat- 
ing the event count. 

For the mass of Galactic bulge stars, a probability density 
in 1k(M/M|^ ) is given by the normalized mass function (e.g. 

lb 



IChabrieili2003h 

*lg(M/MQ)[lg(M/MQ)] 

1.292 exp <! -0.5 



lg{M/MQ)+0.658 



for - 2.0 lg(M/M0) -0.155 .^g) 
0.2546 (M/Mq)-^-^ 

for - 0.155 < lg{M/MQ) 1.8 
for lg{M/MQ) < -2.0 or lg(M/MQ) > 1.8 

which is shown in Figure IA3] 



